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Abstract 

This paper studies false vacuum lumps surrounded by the true vacuum in a real 
scalar field potential in flat spacetime. Fermions reside in the core of the lump, 
which are coupled with the scalar field via Yukawa interaction. Such lumps are stable 
against spherical collapse and deformation from spherical shape based on energetics 
considerations. The fermions inside the lump are treated as a uniform Fermi gas. We 
consider the Fermi gas in both ultrarelativistic and nonrelativistic limits. The mass 
and size of these lumps depend on the scale characterizing the scalar field potential 
as well as the mass density of the fermions. 



1. Introduction 



The spacetime geometry of gravitating lumps is either anti-de Sitter-Schwarzschild 
or de Sitter-Schwarzschild. Dymnikova has discussed the global structure of de Sitter- 
Schwarzschild spacetime, supposing appropriate external matter distribution lj. The 
spacetime geometry of the gravitating fermionic lumps in Ref. [2] is anti-de Sitter- 
Schwarzschild. The mass of such lumps is negative. In this paper we shall investigate 
fermionic lumps with positive mass in flat spacetime. The equations in flat spacetime 
become simple and we expect to obtain similar solutions in curved spacetime where the 
geometry is de Sitter-Schwarzschild. The exploration of the de Sitter-Schwarzschild space- 
time is reserved for future investigation. One may explore such objects in the universe at 
various scales and epochs. 

In field theories with continuous global symmetries, non-dissipative solutions with finite 
energy may exist. These solutions are called non-topological solitons [3J IH|3]. The model 
studied in this paper is very similar to a kind of non-topological soliton called Fermi balls, 
given in Ref. [Hj. In Fermi balls, fermions have zero rest mass and their effective mass is 
due to their Yukawa coupling with a scalar field. In such a model, fermions energetically 
have a tendency to localize at the boundary wall region (domain wall) where their effective 
mass becomes zero. In our model, however, fermions have a nonzero rest mass and their 
effective mass reaches a minimum at the core region of the lump causing the fermions to 
localize at this region. 

In section 2, we set up the problem and introduce useful variables in terms of which 
the field equations are written. In section 3, we investigate the stability of the lumps with 
a positive mass analytically. In section 4, the behavior of the solutions is investigated 
analytically inside the lump, numerically in the boundary wall region, and analytically 
outside the lump. More details of the solutions are given in section 5, with a focus on 
the dependence of the solutions on various parameters of the model. It is seen how the 
solutions change as the energy scale of the model is lowered. In Section 6, we investigate the 
numerical solutions in the transition region with the fermion density having a continuous 
transition to zero. The summary and conclusions are given in Section 7. 
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2. Model with Dirac fermions 



The system that we are considering is a real scalar field coupled with a fermion field 
whose Lagrangian is given by 



We take 



V\ 



£ = iip^d a ip — mipip + -d^ipd^cj) — V[<f>] — g4>ipip ■ 



\{<t> - f2)U 3 - + - \f 2 {h - 2/i)(0 + h) 



(2-1) 



(2.2) 



This potential has two non-degenerate minima at f\ and For further details on the 
potential refer to Ref. In this paper, we employ the natural unit % = c = 1. 

As we shall see below, the scalar field is approximately constant inside the lump so 
that fermions inside the lump may be treated as a uniform Fermi gas. The effective mass 
of the fermion is m e ff = m + g^inside- Depending on the density p n and effective mass m e fj 
of fermions, the Fermi gas can be either relativistic or nonrelativistic. Given p n and m e ff, 
{ipip) = po inside the lump is determined. Let pf be the Fermi momentum of a degenerate 
Fermi gas. Then 



Pn 

Po 



3tt 2 



{pf{p f + m 2 cS ) 1/2 - m ts log 



PF + iPF + m^) 1 ^ 



(2.3) 



2tt 2 «- - - m eS c 

Note that po ~ Pn if fermions are nonrelativisitic. The energy-momentum tensors T^ u for 
spherically symmetric, static configurations are 

1 



V 



Too 
T n 

T 2 2 = T 33 = — - 
others = , 

where primes represent r derivatives. Here 



+ V[<j>]+E f , 

1 „3 



-V\6] + P 



f > 



(2.4) 



F + ™ 2 eS )(PF + ™ 2 e S ) 1/2 



- (m eff ) 4 log 



Pf + {p 2 f + m \ 



(2.5) 
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and 

Pf = Mp^U ~ m «Wr + m ^ /2 + ^ lo § [ PF + (fF + mls)1/2 ] } (2.6) 

are the energy density and the pressure density of a Fermi gas, respectively. T 00 is the 
energy density of the lump which can be used to calculate the total mass of the lump 

M= [ r 4irr 2 drT 00 . (2.7) 



o 



The equation of motion for the scalar field is 

1 d 



r 2 dr 



(ry) + V'[0]+sps(r) = O (2. 



where ps(r) = (tpip )(r), and V'[4>] is the derivative of V[(f>] with respect to 0. 
In order to solve Eq. (|2.8|) we adopt an approximation 

p s (r) = p 6(R 1 -r) , (2.9) 

and we divide the space into three regions: 

I. The inside region (0 < r < R\ ). 

II. The boundary wall region (R\ < r < R2 ). 

III. The outside region (R 2 < r < 00). 

In region I inside the lump, ps = Po and 0(0) is very close to, but still greater than the 
location f s of the minimum of V[4>\ + g<pp : 0(0) = f s + 5(f)(0) with < 5<fi(0)/\f s \ < 1. 
In the inside region, 0(r) varies little from fs so that the equation of motion for can be 
linearized. In region II, ps = 0. In this region the field varies substantially and the full 
nonlinear equation must be solved numerically. In region III, ps = and ~ In this 
paper we focus on the case in which |/i| ~ f 2 ~ / = + / 2 ), f r = (/2 - \fi\)/f < 1 

and /1 — /sr <C / so that the linearization of Eq. ()2.8|) is valid. In Section 6, we will relax 
the third restriction by using a continuous transition from p$ = Po to ps = in region II. 
Schematic behavior of a solution 0(r) is displayed in Fig. ^ 
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Figure 1: Schematic behavior of </>(r). In the lump solutions \f% — fs\ <C |/i| and R 2 — R\ <C 
R,. 

3. Stability of lumps 

As we mentioned earlier, the scalar field inside the lump (region I) is approximately 
constant. This allows us to treat the fermions inside the lump as a uniform Fermi gas. In 
this region, the Yukawa interaction gcfnptp generates an additional linear interaction gpo4> 
for the scalar field. The total potential V[<p] + gpo4> has a minimum at fs < 0. Let R be the 
radius of the lump. To see how false vacuum lumps with fermions become possible from 
energetics considerations, we suppose that the total fermion number Np = 4iTR 3 p n /3 is 
conserved and therefore fixed. It is supposed that V[/i] = e > but the energy density at 
the minimum V[fs] + gpofs = e becomes negative for the lump solution, fs and e depend 
on p n or R. Inside the lump r < R, the scalar field assumes (p ~ fs whereas outside the 
lump <p = f 2 . The total energy of the lump is approximately given by 

E{R) = {£ f ( Pn ,m cS ) + e(p )} V + Sa + E ftWall (R) , (3.1) 

where a is the surface tension resulting from varying cf) in the boundary wall region. V = 
4ttR 3 /3 and S = 47rR 2 are the volume and the surface area of the lump, respectively, po = 
Po{p n , m es) as related by ([2.3)1 . The energy density of a degenerate Fermi gas Sf(p n ,m eS ) 
is given in (|2.5J) . Ef^ wa \\(R) is the contribution to the energy from fermions localized in the 
boundary wall region. Ef >wa ^(R) has been estimated in Ref. [7j to be about 4 v /7riV|^ all /3v / 5' 
for m = 0, where iV/- wa n is the number of fermions confined in the boundary wall region. 
In this paper, we shall focus on the case iV/ )Wa ii <C Np. 

Using Eq. (|2.5p . the energy density of a degenerate Fermi gas can be found in the 
nonrelativistic limit where pp Tn e s, and in the ultrarelativistic limit where pp 3> m e f[. 
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In the nonrelativistic limit 

S f « m cSPn + —p 5 J\ (3.2) 
where A = 3(37r 2 ) 2 / 3 /10, and in the ultrarelativistic limit 

£ f « Bp\l\ (3.3) 

where B = (9tt) 2 / 3 /4. 

For fixed N F , e(p ) ~ e f° r Pn ~ (as R — > oo). Therefore, -E'(-R) ~ (47re /3)-R 3 
for large in both nonrelativistic and ultrarelativistic cases. For large p n (as R — > 0), 
the total potential is approximated by (A/4)0 4 + gp n (p near the minimum so that e ~ 
— ^X~ 1 ^ 3 (gp n ) i ^ 3 . In the case of a lump with a nonrelativistic Fermi gas the £/-term 
dominates over the e-term for small R and E(R) ~ A'R~ 2 , where A' > 0. On the other 
hand, for a lump with an ultrarelativistic Fermi gas the £/-term can only dominate over 
the e-term for small R if B > |A -1//3 g 4 / 3 or 

^ > 3tt 2 . (3.4) 
A 

If Eq. (J23D is satisfied, then E{R) ~ B'R- 1 with fi' > 0. Therefore, with Eq. Q 
satisfied, the total energy -E'(-R) of a false vacuum lump filled with either ultrarelativistic 
or nonrelativistic Fermi gas has a minimum at a certain radius i?, where R is the size of 
the lump. 

Now let us explore some of the conditions and restrictions on the lump solutions in this 
paper. For the lump solutions to become possible we need to have e < 0. Note that e < 
implies that 

PO > ~~~7 ~ ~ ffr ■ (3.5) 

9J 9 

In the examples of the lump solutions below, fs is close to /i, and it is useful to introduce 
the parameter hs = (/i — fs)/ f 1. po and hs are related by 

2A 

Po fh s . (3.6) 

9 

Equations (|3.5j) and ()3.6|) imply that 

2/l5>/r • (3.7) 
Also in this paper |e| ^> eo- In other words 

|e| ~ 9fpo ■ (3-8) 
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First we consider a lump filled with nonrelativistic Fermi gas. As we mentioned earlier, 
the nonrelativistic approximation is justified if m e g 3> pp, which implies that 

m cS » P r-(^f) 1/3 9 f. (3.9) 

We would like to have a lump with positive energy; £f + e > 0. This leads to 

m cS > M „ gf . (3.10) 

Pn 

If the conditions in Eqs. ()3.8|) and (|3.9jl are satisfied, we may have a false vacuum lump filled 
with nonrelativistic Fermi gas with positive mass. Now let us consider a lump filled with 
ultrarelativistic Fermi gas. The ultrarelativistic approximation is justified if m e g <C pp, 
which implies that 

m eS «pr~(^f) 1/3 gf. (3.11) 

For this case also we would like to have a lump with positive energy; Sf + e > 0. This leads 
to 

npo 27r\h s 5/2 ( v 

m ea < —(gf) 1 ■ 3.12 

l e l ' 9 

If the conditions in Eqs. (J3.11)l and (J3.12)l are satisfied, we may have a false vacuum lump 
filled with ultrarelativistic Fermi gas which has a positive mass. 

Note that if the fermion rest mass m >gf, then the effective mass of the fermion inside 
the lump, which is \m — gf\, is smaller than the effective mass of the fermion outside the 
lump, which is m+gf. This condition is satisfied for both cases of lumps with nonrelativistic 
and ultrarelativistic Fermi gas. This is the reason why fermions energetically favor to reside 
at the core of the lump. The situation will change when the fermion rest mass m -C gf or 
m = 0. These scenarios have been studied in the models suggested in Refs. 0| and jHj. 

The stability of a lump configuration against deformation from spherical shape is ex- 
amined in the following way. In Eq. 1)3.1)1 . the surface tension energy and the energy from 
fermions localized at the boundary region Ef jWa u(R) depend only on the surface area of the 
lump and not on its volume. On the other hand, it is clear from Eqs. 1)3.1)1 . ()3.2)1 . ()3.3jl . 
and (J3.8)) that the volume energy will increase if we decrease the volume of the lump by 
deforming it from spherical shape for both nonrelativistic and ultrarelativistic cases. 
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4. Behavior of the solutions 

The procedure to solve Eq. (|2.8jl using the assumption described in the previous section 
is as follows: 
Region I 

First we determine f$ using 

V'[f s }+Po = . (4.1) 

fs and fi are close to each other compared to the length scale /. Denoting 0(0) by 0o, we 
have < (p - f s < 

The equation for 0(r) can be linearized with 8(f)(r) = 0(r) — f$. In terms of z = ru, 

where cj 2 = V'^/s]. The solution to Eq. f|4.2jl regular at r = is 

50(r)=50(O)^^. (4.3) 

As one can see, the ratio 5(f)(r)/5<j)(0) grows exponentially as e z /2z. At the end of region 
I, S(f)/\fs\ needs to be very small for the linearization to be valid. The deviation from fs 
at the origin, S(p(0), needs to be very small to have an acceptable solution. 
The ratio of S(f)'(r) to 8(j>{r) is given by 

5d>'(r) ( cosh z l\ . . , s 

uj ~ uj , (4.4) 



8(f)(r) \sinhz z 
for z ^> 1. 

Region II 

In region II, the equation must be solved numerically because 0(r) varies significantly. 
Nontrivial lump solutions become possible by fine tuning the value of 8<fi(R\). 

How to find the numerical solutions is summarized as follows. First we choose Ri, 
8<p(Ri), and fs- 8<fi'(Ri) is evaluated from the analytic solution in region I. Using the 
boundary conditions 8<p(Ri) and 8<f)'(Ri), we solve the equation numerically in this region. 
The width of the boundary wall region, w — R2 — Ri, is approximately given by 1/y/Xf 

m 

One example of the lump solutions in region II is displayed in Fig. |21 In this example, 
we choose the input parameters g — 1, A = 1, f/Mp\ = 0.0002, f r = Af / f = 0.0002 and 
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Figure 2: <f>(r) of a solution with f/M Pl = 0.0002, f r = 0.0002, A = g = 1, h s = 
(/i ~~ fs)/f — 0.005 and R\/lp\ = 8 x 10 7 . and r are in the units of M P1 and Zpi, 
respectively. 

hs = (/i — fs)l f = 0.005. The resulting output parameters are: po^pi = 8.0593 x 10~ 14 , 
e/M^ = -1.5944 x 10" 17 , |e|/e = 74.738 and cu/M Pl = 2.8495 x 10~ 4 . We note that 
p /p = 0.0101 and e// 4 = -0.00997. For R x /l n = 8 x 10 7 (R^j = 2.2796 x 10 4 ), we find 
a solution with 5(p(Ri)/Mp\ = 5.1491304- • • x 10 -6 . In this example, the value of 8(f) at 
the origin (r = 0) is found from Eq. ()4.3|) to be 1.562 x 10~ 9901 , which explains why one 
cannot numerically integrate 0(r) starting from r = to Ri. A small discontinuity in <fr" 
appears at r = Ri due to the discontinuous change in ps(r). 

Suppose that the fermion mass is 10~ 3 Mpi. Then for the given example above, we have 
a false vacuum lump filled with nonrelativistic Fermi gas which has a total mass of 

M « |e(p n ) + m cSPn + \ n p 5 J 3 } ^— = 1.4 x 10 8 M P1 . (4.5) 
L 10m e fr ' 3 

For the same given example, let us take the fermion mass to be 10~ 19 Mpi, which is approx- 
imately equal to proton mass. In this case, we would have a false vacuum lump filled with 
ultrarelativistic Fermi gas with a total mass of 

M « {e(p ) + pT) ^ = lix 10 37 M P , (4.6) 

In both of these examples, the energy contribution from surface tension a is small (~ 
10 6 Mpi) and can be ignored. 

The behavior of energy density T Q0 is displayed in Fig. El The energy density has one 
sharp peak associated with the rapid variation of <fi. 
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Figure 3: The energy density T 00 for a solution with //Mpi = 0.0002, f r = 0.0002, A = g = 
1, h s = (/i — fs)/ f — 0.005 and R\/lpi = 8 x 10 7 . T 00 and r are in the units of Mpj and 

3(37r 2 ^ 2 / 3 



Zpi, respectively. For r < R±, T 00 ~ e + mp n + 



10m 



p5/3 >0 . 



The field 0(r) approaches /2 for r > R 2 . In the numerical integration, Ri and /g are 
kept fixed while 5<p(Ri) is varied. If 5<p{R\) is chosen to be slightly smaller, <p comes back 
toward, but cannot reach fi, eventually oscillating around = as r increases. If 5<p(Ri) 
is taken to be slightly bigger, then 0(r) overshoots f'2, heading for +00 as r increases. 
With just the right value, 0(r) approaches to f'2 as r — > 00 in region III. Shell solutions 
can appear in which goes back to f\ at r = 00. 

Region III 

In region III, 0(r) is very close to fi- The equation for <ft(r) can be linearized with 
5(j)(r) = f2 — <p{r). The linearized equation is the same as Eq. (J4.2|) with z = ra, where 
a 2 = V"[/2]. The solution in which 5<f)(r — > 00) = is 

e ~Z + Z 2 

5<P(r) = 5<P(R 2 ) ■ —j—. (4.7) 

Z/ z 2 

The value of 8<f>(R2) is determined numerically. 



5. Numerical Analysis in the Nonlinear Regime 

In the previous section, we solved the field equation in three different regions of space. 
In this section, we analyze our numerical results in more detail. 

The lump solutions depend on several dimensionless parameters of our model such as 
f/Mpi, f r , Riu, and pof 3 . The dependence on g and A can be absorbed by rescaling. 
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Figure 4: The / dependence of 5<f>(Ri). A = g = 1, f r = 0.0002, R x u = 2 x 10 4 and 
^5 — (fi — fs)l f — 0.005 are fixed. / is in the units of M P1 . 
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Figure 5: The f r dependence of <50(i?i). X = g = 1, //M P1 = 0.0002, R x uj = 2 x 10 4 and 
= (A - /s)// = 0.005 are fixed. 

When the values of these parameters are varied, the size of the resulting lump solutions 
also varies. In the numerical evaluation, we have taken / and f r in the range between 10 -4 
and 10~ 3 . In addition to these parameters, the fermion mass m e ^ = m — <7|/g| inside the 
lump affects the total mass (energy) of the lump. 

In the numerical investigation g, A, /, f r , R\uo and hs are given to find a desired 
value of 8(j)(R\) for a solution. In Fig. HI 8<j)(R\)/\fs\ is plotted as a function of /. From 
the figure it is clear that the value of the ratio 5(p/\fs\ at R± is independent of /. In 
Fig. 6<f>(Ri)/\fs\ is plotted as a function of f r . 8(f>(Ri)/\fs\ increases as f r increases. 
The numerical evaluation becomes unreliable when 8(j)(Ri)/\fg\ becomes large and the 
linearization in Eq. f)2.8j) becomes invalid. In Fig. 5<fi(Ri) /\f$\ is plotted versus R\uj. 
As R\ui increases, 5<fi(Ri)/\fs\ decreases and approaches a constant (~ 0.016). In Fig. 
6<p(Ri)/\fs\ is plotted versus hs- In this figure, 6<p(Ri)/\fs\ decreases as hs increases and 
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Figure 6: The Ri dependence of 5(j)(Ri) 
h s = (/i-/s)// = 0.005 are fixed. 



A = g = 1, //Mpi = 0.0002, f r = 0.0002 and 
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Figure 7: The h s dependence of S^Rj). X = g = 1, f/M Pl = 0.0002, f r = 0.0002 and 
Rilu = 2 x 10 4 are fixed. 

then begins to increase when hs becomes greater than approximately 0.015. The reason 
for the increase in the value of 5(f>(Ri)/\fs\ is that at some point fi — fs becomes greater 
than 5<f)(Ri), which means that <fi(Ri) < fx- In this situation, the solution will diverge to 
negative infinity if we do not take a large enough value for 8(j){R\). 

In Figs. |H1 El andHUl we plot the mass of the lump M as a function of /, m eff , and 
Riu, respectively. For the data represented in these figures, the nonrelativistic approach 
is justified, i.e. pp m e ff. Combining Eqs. (|3.6j) . (|3.8|) . (|4.5|) . and the fact that 

cu 2 ~2A/ 2 , (5.1) 

we may find a rough estimate of the mass of the lump 

M . ( _ s/+moff) (_^)« ( , 2) 
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Figure 9: M versus m eS . //M P1 = 0.0002, f r = 0.0002, h s = 0.005 and R x u = 2 x 10 4 are 
fixed. M and m e g are in the units of M-p\. 

This rough estimate explains the behavior of M in Figs. |H1 El and Similarly, using 
Eqs. (jSSI), (GEEK, (SSI), (£□]), and the fact that p n ~ 7r(2p /m cff ) 2 / 3 /3 for an ultrarela- 
tivistic Fermi gas, it is possible to find a rough estimate of the mass of a lump filled with 
ultrarelativistic Fermi gas 

M~[-gf+ h s f — 1= \ ; • (5.3) 



^eff / \gV2\J 3 

Using this equation, it is easy to extract the behavior of the lump mass M as a function 
of different parameters of interest. 



6. Continuous transition in the boundary wall region 
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Figure 10: M versus R x u. f/M P} = 0.0002, f r = 0.0002, h s = 0.005 and m efr /M P1 = 10~ 3 
are fixed. M is in the units of Mp\. 
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Figure 11: p s (r) for a/M Pl = 0.01, a/M m = 0.005 and a/M m = 0.001. R x /l vl = 8 x 10 7 
and po/^pi = 8.0593 x 10~ 14 are fixed, ps is in the units of Mp V 

So far we have studied the lump solutions using Eq. (|2.9|) to express ps(r). In this 
approximation, a small discontinuity appears in 0" at r = R x due to the discontinuous 
change in ps at this radius. In this section, we investigate the case in which p$ makes 
a continuous transition from p to in the boundary wall region (region II). To assure a 
smooth transition in this region, we express ps(r) with an equation of the form 

ps(r) = Po { 1 - — \~~T5 TTHTr ■ ( 61 ) 

[ 1 + exp [a(R 1 - r) + 6\ J 

This equation becomes identical to Eq. (|2.9j) in the limit where a — ► oo. In Fig. El we 
plot ps(r) for three different values of a. As a gets smaller, ps make a slower transition 
from po to zero. In Fig. ^1 we show how 5<fi(Ri) decreases as we decrease the value of a. 
In the examples shown in Fig. EI the transition region from p$ = Po to ps = is smaller 
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Figure 12: 8<t>{R x )l\f s \ versus a. f/M Pl = 0.0002, f r = 0.0002, h s = (fi - f s )/f = 0.005 
and R\uj = 2 x 10 4 are fixed, a is in the units of Mpi. 

than the width w = R2 — R\. As a — > 00, the value of in Fig. [T21 converges to 

the value of 50(-Ri)/|/s|, which we find using the approximation in Eq. (J2.9|) . 

7. Summary and conclusions 

In this paper, we have shown that false vacuum lumps with fermions in flat spacetime 
may exist. This configuration is an interesting configuration which may appear in different 
scales and masses. The size of these lumps depends on the typical energy scale of the 
scalar field potential. If the energy scale is very high, the size of the lump is small. If 
the energy scale is low, the lump may extend at a cosmic scale. On the other hand, the 
mass of these lumps depends both on the energy scale of the scalar potential and the 
mass density of fermions inside the lump. For lumps with large masses and sizes, the 
gravitational effects need to be included. The spacetime structure for such lumps will be 
de Sitter- Schwarzschild. We leave the investigation on de Sitter-Schwrzschild gravitating 
false vacuum lumps with fermions to a future project. The stability of the solutions requires 
further examination as well. 

False vacuum lumps with fermions are a type of non-topological soliton that may be 
produced in the early universe due to the spontaneous breaking of an approximate discrete 
symmetry that occurs in a large family of field theories in four spacetime dimensions. It 
has also been shown in Ref. ^9J that the false vacuum lumps of the Higgs field may emerge 
in the early universe in the standard Einstein- Weinberg-Salam theory of electroweak and 
gravitational interactions. In this model, the fermions are the quarks and leptons that are 
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coupled with the Higgs field via Yukawa interaction. It would be interesting to study the 
cosmological effects and applications of such lumps. 
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